Semi-classical calculations of the 
two-point correlation form factor for 
diffractive systems 

E. Bogomolny and O. Giraud 
Laboratoire de Physique Theorique et Modeles Statistiques * 
Universite de Paris XI, Bat. 100 
91405 Orsay Cedex, France 

February 8, 2008 



Abstract 

The computation of the two-point correlation form factor K{t) is 
performed for a rectangular billiard with a small size impurity inside 
for both periodic or Dirichlet boundary conditions. It is demonstrated 
that all terms of perturbation expansion of this form factor in pow- 
ers of r can be computed directly from a semiclassical trace formula. 
The main part of the calculation is the summation of non-diagonal 
terms in the cross product of classical orbits. When the diffraction 
coefficient for the scattering on the impurity is a constant our results 
coincide with the expansion of exact expressions obtained in Ref. Q 
by a different method. 
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1 Introduction 



It is widely accepted that statistical properties of the energy levels of a quan- 
tum problem are determined mainly by the nature of its classical motion. 
Thus spectral statistics of classically chaotic systems follow the random ma- 
trix predictions , [H and spectral statistics of classically integrable models 
agree with the Poisson statistics A lot of numerical evidence (see e.g. 0) 
and partial analytical arguments 0-0 strongly support these conjectures. 

But there are models which (for various reasons) do not fit into the 
above scheme. A noticeable example is given by the so-called diffractive 
systems where one or a few small-size scatters are added to an otherwise 
smooth model. In classical mechanics point-like singularities affect only a 
zero-measure set of trajectories which directly hit them but in quantum me- 
chanics singularities often lead to important effects. Though the addition 
of a scatter to a chaotic system has no effect on its spectral statistics []T0 | 



even one small-size impurity changes completely the statistical properties of 



integrable models ||TT 



The difference between chaotic and integrable models can intuitively be 
understood as follows. In chaotic systems classical trajectories and quantum 
wave functions cover the whole surface of constant energy and the addition of 
a few centers of scattering can not change this property. On the other hand, 
in integrable models trajectories and wave functions are confined to invari- 
ant tori and small-size scatters will, in general, induce transitions between 
different tori, qualitatively changing the nature of motion. 

The investigation of integrable models with short-range scatters has been 
initiated in Ref. |ll| (see also [|12| and references therein) and recently at- 
tracts wide attention from both experimental and theoretical points of view. 
Experimentally, non-hydrogenic atoms in external fields [|1^ are just a possi- 
ble realization of such models. Theoretically, these models are rare examples 
where spectral statistics can be computed analytically. 

In Ref. [0 the two-point correlation function for a rectangular billiard with 
a delta-function potential has been computed exactly and in [|14| the nearest- 
neighbor distributions for this model were obtained. The main conclusion of 
these works is that the spectral statistics of diffractive systems belong to the 
so-called intermediate type. Intermediate spectral statistics describes level 
repulsion like random matrix ensembles but the nearest-neighbor distribu- 
tion has exponential decrease like Poisson statistics. These mixed properties 
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of diffractive systems are similar to spectral statistics of the 3-dimensional 
Anderson model at the metal-insulator transition point \W[ . 



The purpose of this paper is the investigation of spectral statistics of 
diffractive systems by semiclassical methods. Our starting point is the trace 
formula for diffractive systems derived in and which allows to express 
the Green function and all related quantities as a sum over periodic and 
diffractive orbits. The two-point correlation function is obtained as usual by 
smoothing the product of 2 densities of states over a small energy interval. 
It has been shown in [|l^ that the linear on r term in the expansion of 
the two-point correlation form factor K{t) can be obtained by a diagonal 
approximation, which consists of taking into account only contributions from 
pairs of orbits with exactly the same length, and summing them using the 



Hannay-Ozorio de Almeida sum rule [19 



The main difficulty of obtaining the higher order terms of the expansion 
of K{t) in powers of r is the necessity to take into account non-diagonal 
contributions of classical orbits with slightly different actions. The method 
which permits to solve this problem for integrable systems has been proposed 



m 



20|] . Using and simplifying it we construct all terms of perturbation series 
expansion of the two-point correlation form factor K{t) into series of r for 
rectangular billiard with a delta-function impurity for both periodic and 
Dirichlet boundary conditions. The results agree with the expansion of exact 
formulas obtained in 

The model with a delta-function potential is characterized by a constant 
diffraction coefficient. It is of interest to consider more general models with 
a non-constant diffraction coefficient. Though the exact solution of such 
models are not known we construct the perturbation series expansion of 
K{t) in this case as well. 

The plan of the paper is the following. The trace formula for diffractive 
systems is shortly discussed in Section |[ In Section || the detailed discus- 
sion of perturbation series for the form factor of a rectangular billiard with 
delta-function potential and Dirichlet boundary conditions is presented. The 
main part of this section is devoted to the summation of non-diagonal contri- 
butions. In Section ^ the case of Dirichlet boundary conditions is discussed 
and in Section ^ the case of arbitrary diffraction coefficient is considered. 
The expansion of the exact expression of K{t) for rectangular billiards with 
constant diffraction coefficient obtained in is derived in Appendix. 
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2 Trace formula for diffractive systems 



The starting point of the modern semiclassical approximation for multi-di- 
mensional quantum systems is the semiclassical approximation for the (ad- 
vanced) Green function as a sum over classical trajectories with energy E 
starting from initial point x with momentum in the direction n and ending 
at final point y with the momentum in the direction n' |^T|, ||2^ 



Gix, y) = ^ G{{x, n), {y, n')), 

tr 

where the contribution from each trajectory has the form 
G((f, n), {y, n')) = Atr exp{-Sci - i-u). 
Sd = Sci{E,x,y) is the classical action computed along the trajectory, 



(1) 



(2) 



A 



m 



tr 



1 d^Sd 
kk' dydy' 



(3) 



where m is the mass, / is the dimension of the space, y and y' are coordinates 
perpendicular to the trajectory (respectively at initial and final point), k and 
k' are initial and final momenta, and i/ is a phase (the Maslov index) which, 
roughly speaking, counts points where simple semiclassical approximation 
breaks down. Of course, variables (x, n) and {y, n') are not all independent, 
and the sum (|l]) over classical trajectories must take into account, for each 
given starting point and momentum (x, n) , all classically allowed final points 
and momenta. 

For 2-dimensional free motion these formulas take especially simple form 
and the free Green function (in the units m = }^ and h = 1) reads (see e.g. 



tr 



VSnkl 



(4) 



where / is the length of the trajectory. 
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The knowledge of the Green function permits to compute other quantum 
quantities. In particular the quantum density of states 



d{E) = Y,KE - En) (5) 

n 

may formally be written by the means of the advanced Green function as 

d{E) = -^Im j dx G(f , x). (6) 

The contribution from very short trajectories gives the mean level density, 
J, and the integration of (||) over the whole space selects periodic orbit con- 



tributions and determines an oscillating part of level density 

d{E)=d + d^°''\E). (7) 
For any 2-dimensional billiard 

where A is the area of the billiard. The explicit form of d'^°^'^^ (E) depends of 



the system considered (see e.g. pT|-|l23|). 

Diffractive systems discussed in the paper are characterized by the exis- 
tence of singularities which make the classical motion undetermined. Each 
time a classical trajectory hits a singularity there is no unique way to con- 
tinue it. Quantum mechanics smoothes out these singularities and associates 
with each (not too strong) singularity a diffraction coefficient, T>{n, n'), which 
determines a scattering amplitude on this singularity from the initial direc- 
tion n to the final direction n'. In the presence of a singularity the Green 
function in the whole space has two contributions 

G{x, x) = Go{x, x) + Gd(x, x') (9) 

where Gq{x, x') is the free Green function (^) and Gd(x, x') is the contribution 
of trajectories that hit the singularity 

Gd{x, x) = ^ Go{x, (fo, n))V{n, n)Go{{xo, n),x') (10) 
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where T>{n, n') defined by Eq. ( [T0| ) is called the diffraction coefficient. Here 
the quantity Go{x,{xQ,n)) is a contribution to the Green function from a 
classical trajectory starting at point x and ending at the singularity Xq with 
momentum in the direction n; GQ{{xQ^n'),x') is a contribution to the Green 
function from a classical trajectory starting at point Xq with momentum in 
the direction n' and ending at point x' . 

For diffractive systems the density of states (^) can be written as the sum 



of three terms [T^, |T6| 



d{E) =d + dp,o.{E) + dd.o.iE), (11) 

where d is the mean level density, dp,o. is the contribution of periodic orbits 
without singularity, and the third term, dd.o.i^)-, is a contribution from all 
classical trajectories starting and ending at the singularity, called diffractive 
orbits 

°° 1 5 

dd.o{E) = Y] ^^"^ Y] G{ni,n[)V{n[, 712)... G{n^, n'^)V{n'^, Hi), 



^ Tim BE 



(12) 



where G{n, n') is now the contribution to the Green function from a classical 
trajectory starting at the singular point with initial momentum in direction 
n and ending at it with final momentum in direction n' . 

For elastic scattering the diffraction coefficient cannot be arbitrary but 
has to obey the optical theorem which is the manifestation of the quantum 
mechanical unitarity. In the most general form the optical theorem is 



P(n, n') - V{n', n) = -— [ V{n, n")V{n', n") do" (13) 

47r J 

where do" is the angle giving the direction of n". 

In particular, a constant diffraction coefficient should have the form 

V = , (14) 

with real A. 
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3 Periodic boundary conditions 



Let us consider a rectangular billiard with sides a and b and with a point- 
like scatter at a point (xo,?/o) inside the rectangle. We always assume that 
the ratio a^/b"^ is a 'good' (diophantine) irrational number. This generic 
requirement is necessary to avoid 'accidental' degeneracies of periodic orbit 
lengths and under this condition it is possible to prove P that the two-point 
correlation function of such a rectangular billiard in the "free" case (that is 
without scatter) agrees with the Poisson statistics. In this section we impose 
periodic boundary conditions. 

Fig. shows a classical periodic orbit in the rectangle. Since opposite 
sides of the rectangle are identified, one can unfold a classical trajectory to a 
straight line on the plane tiled with rectangles; the images of the point-like 
scatter have coordinates {xq + Ma, + Nb}, M, N & Z. 

m Scattering point 



Periodic orbit (M= 1 , N=2) 

The same orbit unfolded 

Diffractive orbit 



Pencil of periodic orbits 

Figure 1: An unfolded trajectory in the rectangle 

The unperturbed energy levels of a rectangular billiard with periodic 
boundary conditions are 

E = i-ny + {^mr (15) 
a b 

for all integers m,n = 0, ±1, ±2, .... 

When mn ^ energy levels have multiplicity 4. To remove this trivial 
degeneracy it is convenient to consider non-degenerate states only. Asymp- 
totically it gives the factor 1/4 in the periodic orbit density of states and 
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to be consistent with Ref. and Equations (||) and {^U^ it is necessary to 
multiply the diffraction coefficient hj g = 4 (see below). 
Therefore in our case the smooth part will be 

iOTT 

where ^ = a6 is the area of the rectangle. The contribution of the periodic 
orbits to the density of states is given by 



27T y/8Mp 



where 



ip = ^{MaY + {my (18) 



is the length of a periodic orbit. In all integrable billiards periodic orbits are 
not isolated but belong to families which cover an area Ap (see Fig. |ip. For 
the rectangular billiard with periodic boundary conditions the area covered 
by any family of periodic orbits is Ap = A. Since dp,o.{E) corresponds to the 
non-degenerate density of states, the summation in (|1^ is performed over 
all periodic orbits of length Ip with M, > (which corresponds to divide 
by 4 the density (O)); this is expressed by the index . 



3.1 DifFractive density of states 

The third term in the density ( PH ) is the "density of diffractive orbits" 
that is the contribution of classical trajectories starting and ending at the 
singularity. The diffraction coefficient for our model is a constant given 



by (n). The Green functions in (12) can be expressed as a sum over all 
diffractive orbits, that is over all vectors 

Id = {Ma, Nb) (19) 

linking two images of the scatter in the plane tiled with rectangles: 

c.j^anxi-^^^. (20) 
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where Id is given by expression identical to ( [T8| ) 



= vW^FTTW- (21) 

Note that when MN ^ there are (7 = 4 diffractive orbits with exactly the 
same lengths corresponding to ±M, ±A^, As the Green functions enter in the 
diffractive trace formula ( ]I2D always multiplied by the diffraction coefficient 
it is convenient to restrict the summation in (|2^) to M, > and to multiply 
the diffraction coefficient hj g = 4. 

Each term in (|12D involving m Green functions gives the following con- 
tribution to the density of states 



'^■"■^ ' Aimnkdk 
(gV)"" d 

Aim'nkdk ^-^^ ^J^-nkh y/8nkl 

111-- - i^m 

li, . . . ,lm are m diffractive orbits vectors and Ij are their lengths. The Maslov 
index in the periodic case is equal to zero. 

The method which permits to treat such sums has been proposed in [pO| 




+ c.c. (22) 



Its main content is the existence in the sums over orbits Ij of saddle-point 
manifolds which correspond to vectors Ij almost parallel to a given vector L 
(see Fig. D. 

th 

When the m power of the Green function, C", is considered we can 
distinguish different saddle-point manifolds . The first corresponds to the 



case when all m vectors h, . . . ,lm are close to a unique fixed vector L. The 
others are given by contributions where the /j, 1 < j < m, are gathered into 
p groups {2 < p < m) with the following properties: 

(i) within each group all vectors are almost parallel to a vector Li 

(ii) the p vectors Lj are quite far from each other (and in particular their 
lengths are different). 

Let us first order vectors Li, . . . , Lp by a certain manner, e.g. according 
to their lengths 

Li < La < . . . < Lp. (23) 
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o 



Ma 



Figure 2: A saddle-point contribution to G"^{L). 



Because it is assumed that all Lj are different it is always possible. 



Then following ||20[ one can write the formal expression 



P=l Lu...,Lp 



where G is given by (1^) and 



5 ^PJ 



_^ _^ gjfcZi— 3i7r/4 gifcZm— 3i7r/4 



1]^ + . . . + mp — m 

>1 



x5 Zi + . . . + - L] 



^1 ) ■ ■ ■ I^JTi 

. . . 5 f L 



(24) 



(25) 



ip-l+l + • • • + "~ 



for any partition mi + . . . + mp = m of m into a sum of p integers It is 
assumed that vectors Ij all have positive components. N{mi, ...,mp) is the 
number of possible permutations of the Ij and will be discussed later. 

The meaning of this representation is the following (for more details see 
[BO]). The amplitudes corresponding to vectors almost parallel (i.e. belonging 
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to the same group) should be summed coherently but the ones with different 
vectors Lj (that is from different groups) are non-coherent. Therefore after 
a smoothing over a small energy window only the square of the amplitudes 
of non-coherent contributions will survive. 

Let us consider at first the case where m = 2 and p = 1 (see Fig. ^ 

ll,l2 

(in this case the symmetry factor A^(2) is equal to 1). We have to sum 




Figure 3: Contributions to the product of two Green functions 

products of 2 Green functions over all vectors h = {mia,nib) and I2 = 
(m2a, 77,26) verifying the condition 

h + h = L. (27) 

It means that we sum over all positive integers mi,ni,m2,r;,2 such that this 
condition is verified. Changing the sum into an integral, we have to integrate 
over nii and (i = 1, 2). In the coordinates {x, y) where L = (0, L) (see Fig. 
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1^) we set h = {x,y), so that I2 = {L — x,y): since dmidrii = dxdy/A, we get 

/gifcii— 3i7r/4 ^ikl2— Sin/A 
dmidni^== ^— — 6{li + h - L) 

Y pL POO gifcii— 3i7r/4 ^ikh—Siw/i 

= — dx dy , , (28) 

(the condition < a; < L is necessary to have li and I2 in the upper right 
quadrant). In the semi-classical limit where /c ^ 00 we use 

and after performing the integrals 

Gl{L) = / ^ (30) 

More generally this method can be iterated to get 

\ m-l ^ gjfcL-3i7r/4 

(m-1)! VMkL- ^^^^ 

The saddle-point manifold in this case is the space of sets of m vectors almost 
parallel to L such that their sum is equal to L (see Fig. |^). This result can 
also be obtained almost without calculations by using Stokes' theorem. 

It may be noticed that the part of the Green function corresponding to 
the sum of contributions G™, X]m=i(5'^)™~^^]"(-^)5 is the perturbation series 
expansion of the Green function and from (^) one gets 

Y,{9Vr-'Gr{L) = ^M^L)G{L) (32) 



771=1 



where p = 1/A is the density of scatters and 



pikL—Sin/A 
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is the usual contribution to the free Green function of a classical trajectory 
of length L. The result (^) is exactly the exponential attenuation one ex- 
pects for the coherent propagation of a particle in a diffractive medium with 
coefficient of diffraction V. 

The terms . . . ,Lp) can be computed (see ||20|) by the same me- 



thod by considering all possible partitions of the m vectors into p groups of 
almost parallel vectors. But in this case a combinatorial factor which counts 
the exact degeneracies has to be taken into account. For a given partition of 
m into a sum of p integers mi, . . . , irip 

mi + . . . + nip = m (34) 



we obtain a saddle-point manifold in the sum (|25|) by choosing mi vectors 
parallel to Li, m2 vectors parallel to L2, and so on. But the numbers mj do 
not fix a diffractive trajectory uniquely. Trajectories built from the same set 
of primitive diffractive orbits connected in a different order have exactly the 
same lengths and numbers m^. 

All such trajectories correspond to one of the permutation of the following 
sequence of p symbols 

1_^2_^ ...p...p (35) 

mi 5712 mp 

where it is assumed that in the sequence there exist mi symbols of type 1, 
1712 symbols of type 2, etc, and finally m^ symbols of type p. Each symbol 
j corresponds to a different vector Lj defining a certain primitive diffractive 
orbit and a permutation of the sequence ( ^5|) describes how a composite 
diffractive orbit is built. 

For example, at Fig. ^ trajectories with m = 4 and p = 2 are de-pictured. 
Fig. ^ corresponds to all partitions with mi = 3 and m2 = 1 

1112 1121 1211 2111 (36) 

and Fig. ^ represents partitions with mi = m2 = 2 

1122 1212 1221 2211 2121 2112. (37) 

At these figures trajectories of types 1 and 2 are represented respectively by 
horizontal and vertical lines. Since the diffraction coefficient is a constant 
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a) 




b) 



Figure 4: Trajectories with m = 4 and p = 2 with exactly the same length, 
a) mi = 3, 7712 = 1; b) nii = 2, m2 = 2. 

the corresponding saddle-point manifolds (consisting of trajectories close to 
these exactly degenerate ones) all give the same contribution to the density 
dd.o. ', therefore the contribution from one saddle point has to be multiplied 
by their number which is the number of partitions of m into p groups with 
rrii elements in the i*^ group 

(mi + . . . + mr,)! , . 

iV mi, ...,mp) = 38 

mil . . . mpl 

These exact multiplicities of diffraction orbits clearly demonstrate the dif- 
ference between scattering on regular and random arrays. For a random 
distribution of scattering points all trajectories as at Fig. ^ will have differ- 
ent lengths but, because in our model all points of diffraction are copies of 
the same scattering point, these trajectories have equal lengths. Finally 

G^(Li,... ,4)= Yi iiif^) il^^\v ^2^) 



7711 + .. ■-\-'mp—7n ^ i=l ' 

m,- >1 



where G{L) is defined in (p^). 



According to Equation (|22|) , the density of diffractive orbits is the deriva- 



tive of the powers of the Green function. To get the leading term in the 
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semiclassical approximation it is necessary to differentiate only the exponent 
and the diffractive density of states takes the form 

oo 

d,.o.{E) = J2 E d^ll.iLu... ,L,)G{L,)...G{L,) + c.c., (40) 

p=l Li<...<Lp 

where Li, . . . ,Lp are summed over the upper-right quadrant and 

m=p mi + ... + mp^m 

mj>l 

Xf-^)™!-! I (J^)mp~i I Ul) 

where 

R{m,, . . . , m,) = (^Pr + ... + rn,)\ _ 

mil . . . TTLpl 

d!f\ (Li, . . . , Lp) is the contribution of all diffractive trajectories (like those 
at Fig. I) built from mi vectors Li, m2 vectors L2, etc. 
The term with p = 1 is 

<^o (^) = - r -A^r = -(exp(^L) - 1). (43) 

m=l 

Since the contribution of periodic orbits (|l^) can be rewritten 

(44) 

one gets, adding all contributions, 

00 

rf(°-)(E) = E E ,^p)G'(i:i)...G(Lp) + c.c., (45) 

p=l Li<...<Lp 



where 



and dp{Li, . .. ,Lp) = d^fl^Li, . . . , Lp) for p > 2 
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3.2 Two-point correlation form factor 

The 2-point correlation function is related with the level density by the usual 
formula 

R2ie) = {d{E + d{E - '-)) , (47) 

where the brackets denote an energy averaging around E over an interval of 
energy much larger than the mean level spacing 1/d and much smaller than 
energy E. 

The two-point correlation form factor is defined as the Fourier transform 
of the connected part of R2{^) '■ 

K{t) = r ^ (S"''\E + -) d'^°"'\E --)) e^'^^'\ (48) 
J-oo d \ 2 2 / 

(the factors are chosen so that r and K{t) are dimensionless) . 

The advantage of the representation (^5]) is that it is a sum over saddle- 
points (each set of vectors (Li, Lp) for p > 1 is a saddle-point) which are 
uncorrelated, so we can apply the usual diagonal approximation to this sum: 
in the cross product (|48| ) only squares of terms corresponding to a given 
saddle-point will give a contribution. 

Using the expansion of momentum \/ E + t ~ \fE + e/ {2\fE) we get 

~ 1 gie(Li+...+Lj,)/(2fc) 

«^W=E(iS)? E M.(i......yP^^^^ + cx.. (49) 

P=l ^ ' Li<...<Lp ^ P 

where dp{Li, . . . , Lp) are defined in (^). 

As all quantities depend only of the lengths of vectors Li and Li ^ oo 
one can substitute the summation over Li by the integration with the density 

poo 7 

p{l) = / dMdN5{l - ^{Maf + {Nhf) = —. (50) 
Jo 

(Remind that we consider vectors in the upper-right quadrant with positive 
M,iV.) 
Because 

E =^ E (51) 

Li<...<Lp Li,...,Lp 
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for any symmetric summand one obtains 



R2(e) 



dLi 



dL„ 



^ p\ Jo IQAk Jn 16Ak 
+ c.c. 



(52) 



The two-point correlation form factor has the similar form 



p=i 



where 



Kp{T) 



Ank dLi 



dLr 



dpi Jo 16Ak Jo 16Ak 
X 5{Li + . . . + Lp - Andkr), 



\dp{Li 



Mpyui, . . . , Lp)\'^ 



and d = A/lGir. 

The term with p = 1 is especially simple (see (^ 

Ki(t) = I exp(— i r)|^. 



(53) 



(54) 



(55) 



Using the optical theorem ( pTSD this expression can be rewritten in the fol- 
lowing form 



Ki{t) = exp( 



W9_ 
16 



-r . 



(56) 



Substituting (^I]) into (^4]) one obtains the contributions from terms with 
p > 2 (we define Xj by Lj = AnkdTi) 



T 



2 /-oo 



P'- Jo 



dri... dTp J2 Yl 



m,n>p mi+... + mp—m 



ran 



xi?(mi, . . . , mp)R{ni, . . . , Up) 



n 



mi+ni-2 



L(m,-l)!(n,-l)! 



5{ti + ...Tp-r) 



(57) 
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where all mi,ni > 1 . The remaining integral has the form 



J(r) 



dri . . . dTp tI 



. . .T, 



m-p+np—2 



+ Tp - r) 



(58) 



and can easily be computed by the Laplace transform 

j.m+n~p—l 



J{r) 



T 



{m + n — p — 1) 



— ll{mr + nr-2)\. 



(59) 



r=l 



With Eqs. (E 



one obtains the contribution to the form factor 



E E —Rimu...,m,)Rin„...,n,) 

n< ^ — ^ ^ — ^ ^ — ' m,n. 



p\ ^ — ' ^ — ' ^ — ' mn 

m,n>p mi+... + mp = m nj+... + np=r7 
m,- > 1 n,- > 1 



n 

r=l 



{rrir + — 2)! 
imr - DMur - 1)! 



;-8^)- 



{m + n — p — 1) 



-T 



m+n—p+1 



(60) 



Adding all factors together we finally get the complete perturbation series 
expansion of the two-point correlation form factor for rectangular billiard 
with a delta-function impurity and periodic boundary conditions 



Kir) 



|X>|V/16 



(61) 



mnp 



-iVg 



m+p 



p=2 m,n>Q 

where we have defined the rational numbers 



n+p 



T 



m+n+p+1 



A. 



•mnp 



[m + p — l)!(n + j9 — 1)! 
{m + n + p — 1)! 



E n 



G 



mi+rii 



_{mi + l)\{ni + 1)!_ 



(62) 



(the sum is taken over all non-negative integers and Uj verifying mi + 
■ ■ ■ + iTLp = m and rii + ■ ■ ■ + Up = n). Remind that in our case (7 = 4. This 
result coincides with the expansion of the exact formula derived in (see 
Appendix). 

When r ^ 0, K{t) ^ 1 as expected for generic integrable billiards WE\. 
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The special case T> = —Ai obtained for A ^ oo in (JT^) corresponds to 
spectral statistics of the star graphs where Eq. ( |6lD for this value of the 
diffraction coefficient yields 

K{r) = e-'^ + J2 S ^™np(-2r)™+"+P+^ (63) 

p=2 ^' m,n>0 



and can be derived by another method [ p8|| . 



4 Dirichlet boundary conditions 

Let us consider now the case where we impose Dirichlet conditions on bound- 
ary of a rectangular billiard which means that the eigenf unctions of the 
quantum problem vanish on the boundary. It is well known that for this 
boundary conditions the classical motion corresponds to the specular reflec- 
tion on the boundary. For example, Fig. ^ shows a periodic orbit in the 
rectangle. One can unfold the trajectory to a straight line on the plane 





■ 




m 




■ 




^ / 
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■ '' 


■ 


\/m/ 4 




■ 



Scatter 

Periodic orbit (M=l, N=2) 
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Figure 5: An unfolded trajectory in the rectangle 
tiled with rectangles; the images of the point-like scatter have coordinates 
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{(eiXo + 2M)a, (e2?/o + 2A^)&}, where as before (xq, yo) are coordinates of the 
scatter, M, N are in Z and ei and €2 take the values ±1. 
The density of states still contains three terms 

diE)=d + dp.o.iE)+dd.o.iE) (64) 

but in this case the sum (|^) over the eigenvalues is not degenerate any more. 
Therefore the smooth part is 

J=^. (65) 

where ^ = a6 is the area of the rectangle. The density of periodic orbits is 
given by an expression similar to (|1] 



where the summation is still performed over all periodic orbits of length Ip 
with M, iV > 0. But here 



Ip = V(2Ma)2 + (2iV6)2 (67) 

and for the rectangular billiard with Dirichlet boundary conditions the mul- 
tiplicity of periodic orbit lengths is g = 4 for MN 7^ 0. 

The density of diffractive orbits is given by the usual formula (0) with 
diffraction coefficients (0). The Green function in ([l2|) can be expressed as 
a sum over all diffractive orbits, that is all vectors linking two images of the 
scatter in the plane tiled with rectangles: 

where the sum is taken over M, iV G Z since there is no degeneracy of the 
energy levels any more; Id is the length of a trajectory going from (xq, yo) to 
an image of the scatter. 

As in the periodic case we could restrict our investigation to M, > 
but now we have to consider the 16 vectors whose lengths are almost equal. 
Those are vectors linking {xo,yo) to (±2Ma ± Xq, ±2Nb ± yo). The lengths 
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of these orbits are slightly different and can be expanded to the first order 
in 1/M, The vector linking (xq, yo) to (2Ma + eiXo, 2Nb + e2yo) (with 
ej = ±1) has a length 



Id = VC^Ma + (ei - l)xo)2 + {2Nb + (ea - l)l/o)' 

~ /p + (ei-l) — ^ + (62-1)— , (69) 

where lp is given by (^). Taking into account that for this trajectory and 
the Dirichlet boundary conditions the phase exp(— |z7rz/rf) = €162 the terms 
corresponding to a given (M, N) in the Green function ( |68|) can be gathered 
together. There are 16 type of such orbits corresponding to ei,e2 = ±1 and 
all possible changes of signs of M and (see Fig. Their sum is 

16 

J2 e'kk-'^u^ = 16e'^'- sin2(0i) (70) 



i=l 



where 



M 

(f)i = 2k—xoa = kxocosip, (71) 

Lp 

N 

02 = 2k— yob = kyo sin Lp (72) 

bp 

and if is the angle between the periodic orbit vector L and the horizontal 
side of the rectangle. For Neumann boundary conditions sin in the above 
expression should be substituted by cos. 

Therefore the difference between billiards with periodic boundary condi- 
tions considered in the previous Section and the Dirichlet case is the factor 
r{ip) in (|8D 

n,n' p V P 

where the sum is taken over orbits with M, > and (7 = 4. Here we have 
set 

r((/?) = 4sinVi(<^)sin^02(¥')- (74) 
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Compared to the problem with periodic boundary conditions, the Green 
function (^) is modified by a product of factors r{(pi), where ipi is the angle 
between the periodic orbit Li and the horizontal side of the rectangle. One 
can also leave the Green function unchanged but substitute the constant 
diffraction coefficient V by an effective coefficient, S{ip), defined by 

S{v)=Vr{^), (75) 

and taken at the saddle-point. The quantities . . . , Lp) defined by Eq. 

( p5D are unchanged, and performing the same steps as in the previous Section 
one obtains the formula for the density of states quite similar to Eqs. (|4l|) , 

(SIP 



= J2 dp{L,, ... , Lp)GiL,) . . . G{Lp) + c.c, (76) 



where 



d,{L) = ^ exp {gVr{^)L/{2tkA)) (77) 
and dp{Li, . . . , Lp) is given by Eq. (^l]) with 

R{m„...,mp) = (gVr^+-+^r, K + ■ ■ ■ + '^p'^- r{v,rK..r{vpr^ (78) 

mi'. . . . VTipl 

The next steps of the computation of the two-point correlation form factor 
are almost the same as above. The only difference with the periodic case is 
the necessity to know the density of periodic orbits, p(/, v?), with fixed length 
and angle. Like in (pOD one gets 

P(^,V') = ^- (79) 
In the semi-classical limit k oo under the assumption that the ratio 



(xoa) / (yob) is an irrational number the factors 0i and 02 in ([7T|) and ( [f2D are 



equivalent to independent random variables uniformly distributed between 
as that 

(80) 



and gTT. It means that 

< / >= lim - / /(0i(v9),</)2(v2))c?v3 = ^ / / f{(t>iA2)d(t>id(j)2 

fc^oo TT Jq TC^ Jq Jq 



22 



for any smooth function /(0i,02)- 

The integration over the angle therefore gives for any integer n, 



vr 

where we have set 



< riipf >= ^ #i#2 r"(0i, 02), (81) 



4sin2 0isin2 02. (82) 



Of course, the introduction of r(0i, ^2) is just a convenient way to perform 
the generahzed diagonal approximation where one takes into account not 
only trajectories with exactly the same lengths but also ones which have 
equal lengths up to the first order expansion as in Eq. ( pUD (see 
The only necessary integral is thus 



The final expression for the 2-point correlation form factor of the rectangular 
billiard with Dirichlet boundary conditions is similar to Eq. (|61|) 

K{t) = /e-PPWi6\ (84) 



- 1 f-igVY^' figV^ ''^'' 



p=2 m,n>0 



where 



_ (m + p- l)!(n + p- 1)! -A- 



{mi + l)\{ni + l)\ 



(85) 



The sum here is taken as before over all non-negative mi and rij verifying 
mi + ■ ■ ■ + nip = m and rii + ■ ■ ■ + Up = n) and the mean value < . . . > 
is evaluated in (|83D . Note that only the value of A^np differs between the 
periodic case (pT| ) and Eq. (Q); and if r is set equal to 1, ( p5[ ) reduces to 
(|62D. In Appendix it is demonstrated that this expansion agrees with the 
exact results obtained in Ref. fill. 
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5 Non-constant diffraction coefficient 



In the previous Sections we have assumed that the diffraction coefficient that 
appears in the expression of the density of diffractive orbits (|l^ is a constant 
(p!4D for a point-hke scatter. In this Section we consider the general case when 
this diffraction coefficient depends on scattering angles. 

Let a trajectory hit the singularity with momentum in the direction n 
and leave it with momentum in the direction n'. The diffraction coefficient 
V{n, n') is a certain function obeying the optical theorem ([TB|). 



5.1 Periodic boundary conditions 

Trajectories with equal length play a very important role in the computa- 
tion of the form factor. In the rectangular billiard with periodic boundary 
conditions the lengths of periodic and diffractive orbits are given by 

/ = ^{MaY + {my (86) 

and orbits with M, N of different signs are degenerated though geometrically 
they are different. When the diffraction coefficient is a constant all such 
trajectories give the same contribution but when it depends on initial and 
final directions their contribution will be different. 

In the rectangular billiard orbits that hit or leave the singularity with any 
of the 4 angles 

= V, -V, -TT + V?, TT - V2 (87) 

all have the same lengths (see Fig. Let us assume that < (/? < 7r/2 
and label these angles respectively by a = 1, 2, 3 and 4. To take into 
account trajectories with equal length it is convenient to label them by 2 
quantities {ip, a) where the angle ip is uniquely related with the length of the 
trajectory and belongs to the upper-right quadrant, < ip < 7t/2, and the 
index a = 1, 2, 3, 4 describes the form of the trajectory and determines which 
of the 4 angles ( pTj ) is the geometrical angle between the trajectory and the 
horizontal. 

Let us denote by Da/sifjip') the diffraction coefficient V{ga{ip), g^^f')) 
corresponding to a scattering process where the trajectory arrives with angle 
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Figure 6: Four trajectories of same length in the rectangle with periodic 
boundary conditions 

opposite to {ip, a) and leaves the scatter with angle (y?', The explicit form 
of the matrix D{ip, ip') is the following 

D{v,y^')= (88) 

/ V{ip,ip') V{ip,-ip') V{ip,-T: + ip') V{ip,7:-ip') \ 

V{-^,^') V{-^,-^') V{-^,-7r + ^') V{-^,7r-<f') 

'D{—T: + Lp,ip') TT + (p, — vr + — vr + (p') tt + tt — (^') 

y V^ir ^ If, (p') VIt: — ip, ~ip') I?(7r — (p, —TT + (p') 'D{'k — ip^-n — Lp') j 

For instance (see Fig. |^) if the trajectory leaves the scatter with an angle — 
(which in our notation is angle 2), it will come back with an angle vr — so it 
hits the singularity with an angle opposite to a = 2. If it leaves again with an 
angle, say (angle 1), then the scattering amplitude between the incoming 
and the outgoing trajectory will be proportional to V{—ip, ip), which is just 
D2i{(p,(p). The advantage of the matrix representation of the diffraction 
coefficient is that different variables (p describe orbits with different lengths. 
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The scattering amplitude is D 21 



Figure 7: A periodic trajectory in the rectangle with periodic boundary 
conditions 



The proliferation of diffractive orbits with the same lengths is taken into 
account automatically by matrix multiplication. 

We now consider a trajectory made of several periodic orbits. Each time 
the trajectory leaves the scatter with angle a, 1 < a < 4, it comes back 
on the scatter with the unique angle corresponding to an additional phase 
TT, then it leaves the scatter again with any angle (3, 1 < (3 < 4. In the 
general formula ([T^) we sum over all initial and final vectors nij-n'j with the 
only condition that the trajectory is closed, which means that the outgo- 
ing angle of the first diffractive trajectory must correspond to the incom- 
ing angle of the last diffractive trajectory (more precisely they must differ 
by tt). Suppose that we consider a multiple diffractive trajectory consist- 
ing of m diffractive orbits beginning and ending at the singularity and such 
that all unfolded diffractive orbits are almost parallel to an angle (f in the 
upper-right quadrant. Such trajectory will therefore contain a coefficient 
-DoiQ2(v^5 '^)Da2a-i{,'^, V^) • • • -Da„ai(v5, Siuce wc have to take into account 
the degeneracy of the lengths and gather together all terms corresponding to 
trajectories of same total length Zi + ■ ■ ■ + lmi we have to sum over all possible 
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incoming or outgoing angles: 



Note that all angles are equal as we consider trajectories with equal lengths. 
Of course, when "D is a constant tr{D^^{ip, ip)) = {4V)"^ for all integer m > 1. 

When the degeneracy has been taken into account the computation is 
the same as in Section ^ but (gT))"^ is replaced by tr{D"^{ip, </?)) for m > 1. 
The summation of such non-diagonal terms leads to the multiplication of the 
Green function contribution G{L) (given by (|33D ) by an attenuation factor 
similar to the one in Eq. ( |52| ) 



where as above p = 1/^ is the density of diffraction points. 

The higher-order terms in the density of states come from the contri- 
butions of higher-order saddle-point manifolds. Let us consider a composite 
trajectory consisting in m diffractive orbits, among which mi are almost par- 
allel to the direction ipi, m2 to the direction ip2, and more generally are 
almost parallel to the direction (pi for 1 < i < p. Each of such trajectories is 
a permutation of sequence (^) but contrary to the previous Section for each 
diffractive orbit there is an additional label a, = 1,2,3,4 corresponding to 
one of the 4 angles (pTj). 

The total number of different saddle-point trajectories for p, m fixed and 
given nii and angle {ipi, ai), 1 < i < p, is 



the coefficient 4 taking into account the fact that each orbit is 4 times de- 
generate (we will always consider < < 7r/2). 

For instance at Fig. ^ there are 2 families of vectors, parallel to = or 
its images, or parallel to <y92 = 7r/4 or its images. At this figure two diffractive 
orbits are represented: orbits 1,4,5,7,8 are of type 1 and orbits 2,3,6 are 
of type 2. But since now the scattering amplitude T> is not a constant, the 
saddle-points will give different contributions depending on the geometrical 
form of diffractive orbits. 



l<aj<4 




(90) 
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1 8 




Figure 8: Diffractivc trajectories for m = 8,p = 2,mi = 5,m2 = 3 with the 
same length. 

The k-th orbit is determined by a pair {(fii^, ctfe) with ctfe = 1, . . . ,4 and 

For a trajectory with p possible angles, (fi^, can take p different values 
{ifii, . . . , ipp}. Each diffraction on the scatter gives a coefficient which generic 
form is D^j^^ak+ii^ik^ ^ik+i) ^ since the coefficient Dij{ip,ip') is the scattering 
amplitude corresponding to a diffraction where an orbit leaves the scatter 
with the image of type i of the angle (/?, and comes back and leaves again 
with the image j of the angle ip'. In order to take into account all trajectories 
we have to sum over all a^, 1 < < 4, and all ifc, 1 < ^a; < P, for 1 < /c < m 
(of course due to the fact that all orbits are periodic, the index im+i — h)- 

Thus for fixed rrii and (pi the sum over all possible scatterings of the 
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scattering amplitudes gives 
R{{mi, ipi), ... , {nip, ifp)) = ^ tr{D{ipi^, ip,^)D{ipi„ ifi^) . . . D{ipi^,ipi^)) 

(92) 



perm 



where the summation is performed over all permutations of the 

sequence (|35|). 

E.g. for the considered example m = 4 and mi = 1712 = 2 (see Fig. |b 
and (0)) 

-R((2, V^i), (2, (P2)) = tr{D{(pi, (pi)D{ipi, (p2)D{(p2, (p2)D{(p2, (fi)) 
+tr{D{(fi, ip2)D{if2, ^i)D{ipi, lp2)D{lp2, <fi)) 

+tr{D{ipi, ip2)D{(p2, ip2)D{ip2, (pi)D{ipi, (pi)) (93) 

+tr{D{ip2, ip2)D{if2, ^l)D{(fi, (fi)D{Lpi, Lp2)) 

+tr{D{ip2, ipi)D{ipi, ip2)D{ip2, ipi)D{ipi, ip2)) 
+tr{D{Lp2, ipi)D{ipi, ipi)D{ipi, ip2)D{ip2, ip2))- 

The quantity R{{mi, ipi), . . . , (m^, ipp)) replaces the corresponding quantity 
(^2|) in Eq. (^) and reduces to it if "D is a constant. 

The rest of the calculations is exactly the same as above and finally one 
gets the expression (H^) where 



di{L) 



2^ 



(94) 



and dp{Li, . . . ,Lp) is the same as (|4l|) with R{mi, ...,mp) replaced by R 
given by (||). 
Using Eq. ( 



]) and replacing (|^) by (92) we obtain that the two-point 

correlation form factor for the considered case is given by the following for- 
mulas 



Kir) 



--< 



1 + tr ^exp 

00 ^ 



> 



p. 

p=2 m,n>p 



mnp 



m+n—p+1 



(95) 
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with 



A 



;-z/8)™(V8)" 



mnp 



mn{m + n — p — \) 
X < i?((mi, ifi 



w E E n 



(m^ + n— 2)! 
K-l)!(n,-l)! 



(96) 



(mp, ¥?p))i?((ni, (fi),... , {up, (fip)) >, 



where i?((mi, yji), . . . ,{mp,{pp)) is defined in Eq. (^). The sum is taken 
over all partitions of m and n into sums of integers and greater than 
1 verifying mi + ■ ■ ■ + rrip = m and Ui + ■ ■ ■ + Up = n. 

The symbol < . . . > means the integration over all indicated angle vari- 
ables from to it/ 2 



< fi^i,--- ,^k) >-- 



n/2 



n/2 



dipi... I difk f{(pu ■■■ :H^k)- (97) 



5.2 Dirichlet boundary conditions 

Dirichlet boundary conditions give rise to a slight complication as there are 
now 16 diffractive orbits associated with each couple of positive integers 
{M,N) (see Fig. In Eq. (|12D the contribution of a multiple diffractive 
trajectory consisting in m diffractive orbits beginning and ending at the 
singularity and parallel to a fixed direction (f contains products of diffraction 
coefficients V{ipi,ipj) and terms exp{ikli) coming from the Green function 
where li is the length of one of the 16 orbits. If we expand /, for large M, N 
each term exp{ikli) can be written according to ( |5UD as 

^_^^^^ikli-i^Ui _ ^iklp ^i{€i-l)kx(}cosip+i(€2-l)kyosmif {^^) 

where Ip is the length (|67D of the closest periodic orbit and ei,e2 = ±1. As 
in the previous Section we will denote respectively by 1, 2, 3, and 4 the 
angles cp, —ip, — vr + ip and tt — ip corresponding to the symmetry group of 
the rectangle. 

For the periodic case (see Fig. H) a trajectory which leaves the scatter 
with angle [ip, a) comes back to it with an angle which is necessarily the 
angle opposite to the angle a. Therefore for periodic boundary conditions 
the free motion between scatters is fixed and the proliferation of diffractive 
orbits only comes from the scattering process. 
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Figure 9: The 16 trajectories of almost same length in the rectangle with 
Dirichlet boundary conditions. 



For the Dirichlet boundary conditions (see Fig. ^ by small changing of 
the initial angle the trajectory (y), a) can come back to the scatter with any 
of the 4 angles ( |57D corresponding to ip, which leads to another source of 
diffractive orbits with almost the same length. 

Let us denote by Tap^ip) the coefficient corresponding to a diffractive orbit 
leaving the diffractive center with an angle [if, a) and coming back with angle 
opposite to {(p,P). Each Tap{ip) is a product of two factors of the type 



eie2exp [i[ei 



V. 



i{e2 - 1) 



(99) 



where angles = 4>i{(p) were defined in (|7T|) and ([7^) , and = ±1. The 
matrix T = T{if ) describes the free motion modification of the Green function 
and simple calculations show that it has the following form 



/ 1 h ah a \ 

b 1 a ba 

ab a 1 b 

a ab b 1 J 



(100) 
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where a = —e^^'^^ and h = — e^*'^^. 
It can be checked that 



(101) 



where (0i 02 (v^)) is the vector defined by 



V{<P) 



( exp(-i0i - 202) 

-exp(-i0i + i02) 

exp(z0i + Z02) 
\ -exp(20i - 202) 



(102) 



Any diffractive trajectory can be specified by fixing initial and final angles 
for each scattering. If the trajectory leaves the scatter with an angle ((/?, /c) 
and comes back to it with angle (</?, I + tt) (that is, an angle opposite to I) it 
is attributed a coefficient Tki{ip) ( |100| ). If then it scatters to an angle {(p',j) 
it gets the coefficient Dij{ip,ip') defined in (^) and so on (see Fig. |1^). For 




Figure 10: An example of periodic trajectory for m = 3 

example, a diffractive trajectory consisting in m diffractive orbits parallel to 
the direction ip will be associated to a pre-factor 

Taia{Da'^a2Ta2a'2 ■ ■ ■ Ta^^a'^D^'^a-^ = tr(TD)"'. (103) 
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As in Section |^ assuming that the ratio (xqci) / (yob) is an irrational num- 
ber then as — >■ oo angles 0i and 02 act as independent random variables 
uniformly distributed between and 11/2. 

To take into account higher-order contributions we have to sum over all 
possible types of vectors (ctj, a'j = 1 to 4) and all possible angles ipi^ G 
{ipi, . . . , ipp} in terms like 

= ir{D{(p,^,ifi,^)T{if,2)D{!fi2,(pi.^)T{!fi.^) ...) (104) 

under the constraint that each angle (fi appears exactly times. Using 
the property ( |101|) one can rewrite the above expression by introducing the 
effective diffraction coefficient S defined by 

Siif, = ^')Vpi^') (105) 

which describes all possibilities of scattering and free motion between the 
scatter. Since S{(p, ip') is a number we get 

= Siipi,, (^iJ5((^i2, <pi,) . . .S{<pi^, <pi,) (106) 



Let us define a quantity similar to Eq. 

perm 

(107) 

where the sum is taken over all permutations (ii, ^m) of the sequence (p5|) 
and it is implicitly assumed the set ipi contains mi terms of cpi, m2 terms of 
ip2,- ■ ■ , terms of ipp. It can be checked that Eq. (|107| ) reduces to (^) if 
"D is a constant, since we have 

J2%iip)V,iip) = griif) (108) 

with (? = 4 and r((^) given by (|7^) . The final expression for the two-point 
correlation form factor, K{t), follows from exactly the same consideration as 



33 



in previous Sections. It reads 



where 

A. 



K(t) 



--<< 



exp 



-I r 



>> 



oo ^ 



p. 

p=2 m,n>p 



(109) 



mnp 



(-i/8)™(z/8)" 
mn{m + n — p — 1)\ 



E En 



{rrii + ni- 2)! 
fm,- - DUn, - 1] 



X << i?((mi, v^i), • • • , (mp, V3p))i?((r2i, v^i), • • • , (n^, v^p) >> . (110) 

As usual the sum is taken over all partitions of m and n into sums of p 
positive integers. The symbol << ... >> denotes the double average over 
random phases (pi as in Eq. (^) and geometrical angles ipj as in Eq. (pTf). 



6 Conclusion 

We have discussed the spectral statistics of a rectangular billiard with a 
small-size impurity inside and developed the method which permits the ex- 
plicit construction of perturbation expansion of the two-point correlation 
form factor for this system, K{t), into series of r 

Kir)=J2cnT\ (111) 

n 

Using the method of Ref. we demonstrate that after the summation over 
diffractive orbits with approximately the same lengths the oscillating part of 
the density of states can be written in the form ( ^5]) 

oo 

d(°-)(E)=^ d,{U,... ,Lp)G{L^)...G{L,)+c.c., (112) 

p=l Li<...<Lp 
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with 



dp{Li, . . . ,Lp) — — — — - — - y2 — 

X ^ i?((mi,(^i),... ,(mp,(^p)) 

m]^ + ... + mp— m 

TTlj > 1 

^2tkA' im^-iy/"^2tkA' K-1)!' ^ ' 

where the summation is taken over all partitions of m into sums of p posi- 
tive integers m = mi + . . . + nip and i?((mi, ipi), . . . , {nip, ipp)) ( |107| ) is the 
contribution of all diagrams describing the scattering process of composite 
diffractive orbits consisting of mi orbits of type 1, m2 orbits of type 2, and 
so on. 

To construct such diagrams it is convenient first to write down all permu- 
tations of the sequence (^) of m elements with m^ elements of type i. The 
total number of these permutations is 

N{mi,...nip) = — -. (114) 

mi! . . . mpl 

Each symbol i represents the angle ipi and the constructed sequence of sym- 
bols is considered as the representation of the scattering process. The sym- 
metry group of the rectangular billiard leads to the necessity to consider 
together all symmetry partners of a given trajectory, that is all possible 
choices for each diffractive orbit. 

For periodic boundary conditions there are 4 different possibilities (^) 
as indicated at Fig. ^, whereas for Dirichlet conditions it gives rise to 16 
trajectories with approximately the same lengths as at Fig. ^ To take into 
account all these degenerate and almost degenerate trajectories it is useful 
to use the matrix formalism developed in the previous Sections. All the 
scattering processes between 2 consecutive angles ip and ip' are described by 
the matrix D{if,if') defined in Eq. (^) for periodic boundary conditions, 
and by the quantity S{ip, ip') defined in Eq. ( |105|) for Dirichlet ones. 

The quantity i?((mi, (pi), . . . , {nip, ipp)) is the sum of such contributions 



for all permutations of our standard sequence (|35|), that is all possible order- 
ings of the sequence of diffractive orbits. 
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The knowledge of i?((mi, ipi), . . . , (m^, ipp)) permits easily the computa- 
tion of the two-point correlation form factor in all cases. The advantage of 
the representation ( |112D is that different terms in it are non-coherent and the 
two-point correlation function is simply given by the sum of the squares of 
each terms as in Eq. (^). This 'generalized diagonal approximation' leads to 
explicit formulas for the form factor Eqs. (|95D, (|96|) for the periodic boundary 
conditions and Eqs. ( |109| ) and ( |110| ) for the Dirichlet ones which solve the 
problem of the summation over non-diagonal terms for integrable systems. 



Appendix 

The exact expression for the two-point correlation function, R2{^), for a rect- 
angular billiard perturbed by a small-size impurity with constant diffraction 
coefficient has been derived in Ref. The aim of this Appendix is to obtain 
its perturbation expansion for large e and the corresponding expansion of the 
two-point correlation form factor, K{t), into series of r. 

According to Eq. (143) in |]l| (with slight changing of the notations) the 
two-point correlation function can be written as multiple integrals 



oo rco 



i?2(e) - 1 = / da, / da2 e-"^<''(™i'™2)> (115) 
Jo Jo 

x[< rJo(2rv/SIai)e*"("i+"2) >^ + < r Ji(2rySI^)e^"("i+"2) + c.c. 
where 

■ ^oo 



J(ai, aa) = (ai - a2)(l + — : + 2ie"^^ / M2^/a2t)e'^dt) 

nv 



- z[«iJo(2v/^) + v/^Ji(2v/^)]e*("i+"^) (116) 

Here Jo{x) and Ji{x) are the usual Bessel functions. The renormalized cou- 
pling constant v' is connected with the diffraction coefficient, T>, by the re- 
lation 

4 Anv' 
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The variable r = r(0i, ^2) depends on the boundary conditions (see 

{1 for periodic boundary conditions 

4sin^((/)i) sin^(02) for Dirichlet boundary conditions , 
4cos^((/)i) cos^((/)2) for Neumann boundary conditions 

(118) 

where 0i and 02 are independent angles distributed uniformly between and 
71/2. The symbol < . . . > denotes the mean values over these angles 

< /(01, 02 >= \ rf01 ^^02/(01, 02). (119) 

Jo Jo 

In particular we have < r(0, 0) >= 1. Changing variables e to x and ai, a2 
to ai/x,a2/x, one gets 

poo pcx> 

i?2(x) - 1 = / t^ai / rfa2e-"^<^^^'^)> (120) 
Jo Jo 

X < r Jo( — )e - >^ + <rJi( — )e - + c.c. 

a; X 

To compute the expansion of the above expressions into powers of 1/a; it is 
convenient to use the transformations proposed in |p7|. 
From ( 116 ) it follows that 

+ ^\xJ[—,—) = 2Jq[-^ )e - (121) 



dai da2 I x x x 



and 



^xJ(^, ^) = -2.^J,(^^)e^^. (122) 

ox X X X X 

Using the fact that Jq{x) = —Ji{x) and J[{x) = Jq{x) — Ji{x)/x, one can 
check that 

^ +^VxJ(^,^) = ^|^xJ(^,^). (123) 



dai da2 J x x a\a2 dx^ x x 
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Therefore, 



d d 

+ 



aia2 dx"^ 

d d 



,7/ roL-i rctn ^ 



TT 



(9a;i da2 



rai ra2. 
X < J{ , ) > 



x^ / d ^,rai ra2s 



X X 

2" 



(124) 



^_^a:<J(-^,-^)> 



and according to Eqs. ( |121| ) and ( |122|) the previous equation can be transform 
as follows 



47]-2 



d d 

+ 



x^ (9^ 



dai da2 J a\a2 dx'^ 



(125) 



< rJo(C)e''^"^ + < rJi(C)e*'^"^ >^ 



g-.x<j(ifi,i^)> 



where C = 2r^/aia2/x. 

The right-hand side of this expression is the pre-factor which appears 
in Eq. ( p.20|) and, consequently, one can express the two-point correlation 
function in the form 



R,{x) - 1 



daida2 
47r^x^ 



d d 
+ 



x^ (9^ 



dai da2 J (X\a2 dx^ 



-|- c.c. 



_^a:<J(— 1,— ^)> 



(126) 



The derivatives with respect to cti and 02 can be computed by integration 
by parts. It gives 

d f d d 



^ dai \dai ' da2 



+ ^ I e— <j(i^,i^)> 



27r(re*™^/^)e~™^(^"^) (127) 



and 



da — (-^ 
^ (?a2 ^"^2 



_,rx<J(— L,— a)> 



2vr(re^""i/^')e-™^(^+^). (128) 
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Therefore 

^Q^i^Q^^ f_d_ _d_\ ' 







^^(^ + ^^^) ^+c.c. (129) 



We are interested in the Fourier transform 

POO 

K{t) = dx {R2{x) - l)e2^--^ (130) 



of the two-point correlation function. 

The term ( |129| ) gives the following contribution 

Ko{t) =< 6^*^™+ + e^^'^™- > -2. (131) 

where 

^± = , / . (132) 

are the poles in (^). Using the definition ([1171) of V, Eqs. (|131|) and ( |13^ ) 
lead to 

irX>q ir'Dq 

Ko{t) =< e"— " + e— " > -2, (133) 

where we introduce g = A to he consistent with notations in Section |^. 

The other contributions come from the second derivatives with respect to 
X in ( 11261 ). Using expansion (112) of Ref. [|I| one gets 



X < J[ , )> = ai + ^2 H :(«i-«2) (134) 

X X TTV' 



oo 



^. (2ai)"+i(ia2)"+^ {m + n)\ < r™+"+2 > 
* j^m+n+i m!(m + l)!n!(n + 1)! ■ 

m,n=0 ^ ' J \ ' J 

The expansion of the exponential of this quantity yields 

p=0 ^' m,n>0 

"•1 H hJip— 
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As 



a 



a+l ' 



(136) 



the integration over ctj gives the expansion of the two-point correlation func- 
tion into series of 1/x 



■p\X) , 



(137) 



where Ro{x) is given by ( |129| ) and 

(2i7r)P ^ {m + n + p + ly.Amnp 



Rp{x) 



pi 



m,n>0 



16tt 



m+p 



Vg_ 

IGvr 



n+p 



where 

A, 



mnp 



[m + p — l)\{n + p ~ 1)! 
{m + n + p ~ 1)\ 



E n 

nii ,nj >0 i=l 



mi+rii 



^ j,mi+ni+2 ^ 



{rrii + ly.im + 1)1 



;i38) 



. (139) 



The expression ( |139| ) for Amnp is the same as in Eq. ([85|). 

Taking the Fourier transform (|13CI|) of ( p.38| ), we get the corresponding 
expansion of the two-point correlation form factor. The term corresponding 
to p = 1 can be transformed to the form 



i^i(r) =< (e-'^^-l) (e 



irT>g 



1 >. 



(140) 



The sum of the contributions (|133|) and ( [140| ) (plus a 1 coming from the 
(5-function in R2{x)) gives 



irVg irT>g I'D] gr 

< e 8 8 ^ >=< e is ^ > . 



The Fourier transform of the sum ( |138|) gives the terms for p > 2 



p=2 m,n>0 



m+p 



n+p 



T 



m+n+p+1 



(141) 



(142) 



which coincides with the result (84). 
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